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ON PERMANENTAL POLYNOMIALS OF CERTAIN 
RANDOM MATRICES 



Abstract. The paper addresses the calculation of correlation func- 
tions of permanental polynomials of matrices with random entries. 
By exploiting a convenient contour integral representation of the 
matrix permanent some explicit results are provided for several 
random matrix ensembles. When compared with the correspond- 
ing formulae for characteristic polynomials, our results show both 
striking similarities and interesting differences. Based on these 
findings, we conjecture the asymptotic forms of the density of per- 
manental roots in the complex plane for Gaussian ensembles as well 
as for the Circular Unitary Ensemble of large matrix dimension. 



The permanent of an N x iV matrix H with entries (i,j = 
1, . . . , N) is defined by 



where the sum is taken over all permutations a of {1, ... , N}. In strong 
contrast to determinants, computing permanents is ^-complete prob- 
lem jSHl, and a considerable effort was spent on developing various 
approximation methods, see [2JJ and references therein. Permanents 
have important combinatorial meaning, in particular the permanents 
of (0, 1) matrices enumerate matchings in bipartite graphs, see e.g.[20j. 
Permanents also have applications in physics of interacting Bose parti- 
cles, see e.g. |51j . A concise introduction into properties of permanents 
can be found in [3H], the standard reference is [32] . 

The permanental polynomial afforded by H is defined as the perma- 
nent of the characteristic matrix, i.e. 



(1.2) p (//) = Per (/il N - H) = fi N -a 1 fi N - 1 +a 2 fi N ~ 2 -. . .+{-l) N a N 
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where Ijv stands for the N x N identity matrix. In particular, it is 
easy to show that a\ = TrH and = Per(H). Although a per- 
manental polynomial is not preserved by similarity, it is preserved by 
permutational similarity: 

(1.3) Per (//ljv — P -1 HP) = Per Guljv — H) 

for any N x N permutation matrices P. Permanental polynomials at- 
tracted some interest in graph theory and its applications to chemistry, 
see e.g. [131 ESI US]- In particular, some general statements on location 
of the roots of permanental polynomials (known as permanental roots) 
can be found in |3H1 ESj • Some curious series expansion of moments of 
permanental polynomials were developed in [TT] . 

The goal of the present paper is to initiate research on permanental 
polynomials of various classes of random matrices. Although some 
questions concerning permanents of random matrices were investigated 
earlier, see the books [33113 as well as more recent results in |4*51 I46j . 
the present author is not aware of any related studies on permanental 
polynomials. The final goal of the whole project should be finding 
various statistical characteristics of the permanental roots, the mean 
density in the complex plane being the simplest, yet nontrivial example. 
Formulated as such the problem seems to be quite challenging, and no 
results are currently available. 

To get some insight into the problem in the present paper we address 
simpler, yet informative objects which we call the correlation functions 
of permanental polynomials. The n-point correlation function depends 
on n complex parameters . . . , \i n and is defined as the expectation 
value 

(1.4) (Per (^In ~ H) Per (fi 2 l N - H) ... Per (fi n l N - U)) H , 

where the angular brackets stand for the averaging over the distribution 
of the matrix entries of H. Particular classes of the distributions we 
are dealing with will be discussed explicitly later on. 

To this end it is appropriate to mention that in the last few years 
it became clear that the general correlation functions of characteristic 
polynomials d(fi) = det (/iljv — H) of random matrices of various kinds 
are extremely informative objects, with a rich mathematical structure 
and important applications. One of the sources of interest in such 
objects originated from attempts by Keating and Snaith[33J EH E2] 
to get an analytical insight into statistical properties of the Riemann 
zeta-function and related functions, whose zeroes are believed to share 
statistically many features with the eigenvalues of unitary (or Hermit- 
ian) random matrices. Brezin and Hikami [TTJ ^] understood that 
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the correlation functions of products of characteristic polynomials for 
invariant ensembles of random matrices can be represented as determi- 
nants made of orthogonal polynomials generated by the random ma- 
trix measure, see (|1.10|) below. Fyodorov and Strahov generalized this 
result to correlation functions involving both products and ratios of 
characteristic polynomials |28| 1571] (an elegant proof is by Baik, Deift 
and Strahov [5|, and important extensions to other symmetry classes 
are due to Borodin and Strahov, jT0|) . The latter objects are especially 
important for applications in physics of quantum chaotic systems, see 
e.g. Andreev and Simons[3], and more recently in Quantum Chromo- 
dynamics , see e.g. Verbaarschot and Splittorff j^T], Fyodorov and 
Akemann J2U 12] and references therein. Characteristic polynomials 
of non-Hermitian/non- unitary random matrices of several types were 
studied most recently by Akemann and collaborators j3J and by Fy- 
odorov and Khoruzhenko |27j . 

All these facts providing an additional strong motivation for the 
research, a number of far-reaching generalizations and extensions of 
the above-mentioned results were obtained in recent years by various 
groups, see [Si HH UH UB UJH EH]- In addition, nice combinatorial inter- 
pretations of moments of the characteristic polynomials were recently 
revealed in H2| E3 H2| • 

A major part of the recent progress in dealing with the characteristic 
polynomials was essentially possible due to the fact that those poly- 
nomials depended only on the eigenvalues Ai,...Ajv of random ma- 
trices, not on their eigenvectors. The joint probability density (j.p.d.) 
V(\\, . . . , Ajv) of those eigenvalues is well-known for the standard classes 
of random matrices. 

As one of the most important examples, consider the space M. of 
N x N Hermitian matrices H = (i/y) = H*, with the probability 
measure on M. being chosen according to 

N N 

(1.5) V{H)dH = ce-^^HdHuHdHSdHi 

i=l i<j 

where Hij = + iHj- in terms of its real and imaginary parts, and 
c stands for the corresponding normalization constant. The function 
V(H) is known as the potential and is usually assumed to be a polyno- 
mial in H of even power with real coefficients, the coefficient in front of 
the highest power being positive. Such potentials define the so-called 
Unitary Ensembles [18J, which is a distinguished subset of distributions 
invariant with respect to unitary conjugations: H — > UHU*, with U 
standing for any N x N unitary matrix from the group U(N). Let us 
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further use the spectral decomposition H = UAU* of the Hermitian 
matrix H in terms of the diagonal matrix A = diag(Xi, . . . , A at) of its 
(real) eigenvalues, and some U G U(N). This decomposition induces 
the corresponding decomposition of the measure (see (HI 123 HI] for 
more detail): 



A? 



1.6) V{H)dH = C N e- NTrViA) A 2 {A) J] d\ dfi N {U) 



i=i 



where dn^iU) stands for the corresponding invariant (Haar's) measure 
on U(N) (normalized to unity), A(A) = Elio^i — ^j) * s the so-called 
Vandermonde factor, and Cat is the corresponding normalization con- 
stant. Integrating out the [/-variables immediately gives j.p.d. of real 
eigenvalues in the form 

N N N 

(1.7) P(Ai, . . • , Aat) d\i = C N HiXi-X^HdwiXi) 

1=1 i<j i=l 

where the measure dw(X) is related to the potential V(X) as dw(X) = 
e~ NV ^dX. It is then a straightforward exercise [T%| to show that 
the expectation value (d(fi))n of the characteristic polynomial is simply 
given by 

(1.8) (det(jil N -H)) H = Tr N (jj) 

where 7Tfc(x) stands for the k— th monic orthogonal polynomial ( tti~(x) = 
x k + lower powers), with the orthogonality being understood with re- 
spect to the measure dw(x) on the real axis: 



(1.9) / 7r j (x)7r k (x)dw(x) = CjC k 5 jk 

In fact, the relation equivalent to ()1.8|) was known in classical theory 
of orthogonal polynomials (although, without any reference to random 
matrices) for more than one hundred years P=8 . It is its modern random 
matrix interpretation that made evident a deep relation between the 
random characteristic polynomials and orthogonal polynomials, and 
enabled those far-reaching generalizations that were already mentioned 
above. In particular, the simplest case of the Brezin-Hikami result reads 



1.10) fe)«^ = — det' * nM MfM2 ' 

/ii — /i 2 
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Naively one may expect that nothing similar is to be valid for per- 
manental polynomials, precisely because of the lack of invariance under 
general similarity transformations. This fact makes seemingly imprac- 
tical presenting permanental polynomials in terms of the matrix eigen- 
values, and forces us to search for an alternative technique of evaluating 
the random matrix averages featuring in (jl.4j) . both over the eigenval- 
ues and the corresponding "angular" variables (=eigenvectors). To this 
end, we start with revealing a possibility to represent permanents as 
a multiple contour integral, see Lemma (|2.1j) . Such a representation 
turned out to be very useful, and allowed us to calculate explicitly the 
expectation value of the permanental polynomial for a general Unitary 
Ensemble of Hermitian random matrices. In fact, integrating out the 
"angular" variables can be done explicitly with the help of a variant 
of the famous HarishChandra-Itzykson-Zuber (HCIZ)formula[30j 152] . 
and is based essentially on interpreting that formula in terms of the 
Schur function expansion method, see e.g.|H|. It came as a some- 
what surprising result that the expectation value of the permanental 
polynomial can be generally expressed as a one-fold integral of the same 
monic orthogonal polynomials featuring in ()1.8)) . see (Jl.llj) . In fact, we 
prove the following 

Theorem 1.1. The expected value of the permanental polynomial for 
random matrices taken from a Unitary Ensemble characterized by a 
potential V(x) is given by 

(1.11) ( p ( fl ))=a N 1 J e- NV ^(fi- A) 2Af " 1 7r Ar _ 1 (A)dA 

,, k r, (2iV - 1)! [ -NV(\) \2N-l-k„ A\ j\ 

= a N 2_^a k fi, a k = — J e 'X 7rjv-i(A) dX 

where 7Tat(A) = X N + lower degrees is the N—th monic polynomial or- 
thogonal with respect to the measure dw(x) = e~ NV ^dx on the real 
axis. 

For a particular case of the Gaussian Unitary Ensemble (GUE) char- 
acterized by the measure V(x) = x 2 /2 the integral in 1)1.11)1 can be 
evaluated explicitly yielding , even more surprisingly, the relation al- 
most identical to (J1.8)) . with the only simple change /i — > — ifi, see 

(1.12) (p(^))gue = i N (d(-ifi)) G UE, d(n) = det (/J,1 N - H) 

Note: In this work we consider general Hermitian matrices with 
the eigenvalues A, G (—00,00). Sometimes, however, one may wish 
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to deal with positive definite Hermitian matrices, and thus to confine 
the eigenvalues to the positive semi-axis Aj G [0, oo) (e.g. the so-called 
Wishart (or Laguerre) Ensemble, see [22], with V(X) = A, A > 0). 
Further choices confining eigenvalues to an interval of the real axis are 
also possible. Extensions of our result to all those cases are self-evident. 

The explicit evaluation of the lowest non-trivial correlation function 
of permanental polynomials turns out to be possible also for one more 
important class of random matrices, those from the Circular Unitary 
Ensemble (CUE), see Section |3] The corresponding matrices are gen- 
eral unitary N x N and are considered to be uniformly distributed 
on the group manifold according to the corresponding Haar's measure. 
The calculation is based on a recent progress in evaluating some inte- 
grals over U(N) [23132], and again has the Schur function expansion 
in its heart. 

Although in principle the Schur function expansion method can be 
extended to obtain higher-order correlation functions, the actual cal- 
culations for general invariant ensembles as well as for the Circular 
Unitary Ensemble become cumbersome and the progress is yet to be 
achieved. We therefore concentrate on two particular cases where a 
progress can be achieved by a different method: the Gaussian Unitary 
Ensemble of Hermitian matrices (Section ^ and the so-called Ginibre 
Ensemble (Section IHJ) of general complex matrices with independent, 
identically distributed normal entries. Our starting expression is the 
same contour integral representation of a permanent, but the specific 
choice of the measure allows to perform the ensemble averaging in a 
closed form without resorting to HCIZ formula. Further progress is 
based on a method frequently employed in the physics of disordered 
systems and known there as the Hubbard- Stratonovich transformation, 
see e.g|21]. In this way we managed to represent the general n— point 
correlation function via n— fold integrals: 

Theorem 1.2. For any integer n = 1,2,... the general n— point cor- 
relation function of permanental polynomials of GUE matrices has the 
following integral representation: 

(1.13) (flpM) = f e"^ Trq2 [Per (M - q)] N dq 

\k=l I GUE J 

In this expression M = diag(fii, . . . , fi n ) is n x n diagonal matrix and 
q is a general n x n Hermitian matrix. 
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Note. For fii = fi 2 = ■ ■ ■ = = ^ the identity ()1.13j) can be 
written as 
(1.14) 

J e -% TrU2 [Per (jjl1 n - H)] n dH = J e~^ Trc ? [Per (jil n - q)] N dq 

Similar identities are known to hold for moments of GUE characteris- 
tic polynomials, and are frequently called "the duality relations", see 
e.g. |2E1 122] and references therein. In that context they were com- 
monly thought to be intimately connected to the underlying integrable 
structures (Toda lattice hierarchies, see e.g. [22!)- The fact that they 
emerge in the present context may indicate some hidden integrability 
lurking behind the problem of permanental polynomials. 

The formula ()1.13j) is already a substantial simplification, especially 
if we have in mind random matrices of large size N, so that n can be 
small in comparison with the size of the matrix. Nevertheless, to get a 
closed form result for any n and ./V amounts to evaluating those inte- 
grals and remains an outstanding problem. However, for the two-point 
correlation functions n = 2 the integral can be evaluated explicitly, and 
the result is expressed in a very attractive "determinantal" form which 
is actually again almost identical to (jl.lOj) for characteristic polynomi- 
als. Namely, we found that 

(1.15) ipM pM)gue = (d(-ifii) d(i^ 2 ))GUE ■ 

We also show that a formula similar to (jl,13|) can be derived by essen- 
tially the same method for the Ginibre Ensemble, see ()6.3|) . Finally, the 
same method works also for the ensemble of real symmetric matrices 
with the Gaussian probability measure (the so-called Gaussian Orthog- 
onal Ensemble, GOE), for which we again obtain explicit expressions 
for the one- and two-point correlation functions of permanental poly- 
nomials, see SectioiQ Close similarities with the corresponding results 
obtained earlier[12j for characteristic polynomials of GOE matrices are 
again apparent. In particular, both the relations (|1.12|) and (jl,15j) 
remain valid. 

The last (but by far not the least) point to be mentioned is that one 
is usually interested in the asymptotic, large-iV behaviour of random 
matrix characteristics. In particular, for GUE characteristic polynomi- 
als asymptotic behaviour is very different for \i e (—2, 2) (the so-called 
"bulk scaling" regime), for = 2 + 0(iV- 2 / 3 ) ("soft edge scaling") and 
for all other complex values of \i JH]- This is naturally related to the 
fact that all N eigenvalues \ of GUE matrices, being real, concentrate 
in the limit iV — > oo in the interval A G [—2, 2] with probability tend- 
ing to one. The corresponding limiting density of eigenvalues is given 
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by the famous Wigner semicircular law: Poo(A) = — A 2 . More- 

over, within that interval those eigenvalues are non-trivially correlated 
on the distances |Aj — Xj\ of the order of 0(1/N). The latter fact is 
reflected, in particular, in the existence of the following limit, [H], see 
also |23: 
(1.16) 

(d(fa) d{fa)) GUE smx 

N^oo {d 2 {fa)) GUE X 

whenever fix, fa £ (~ 2,2), hin^^oc N(fa — fa) < oo. Moreover, the 
limiting expression (frequently called the "Dyson kernel") is not re- 
stricted by GUE, but is valid for a very broad class of Unitary Invariant 
Ensembles of random Hermitian matrices ("universality"). 

Adopting those facts to the case of GUE permanental polynomials 
by the relation ()1.12|) we see a special role of the interval /i G [—2i, 2i\ 
of the imaginary axis. Combining ()1.15|) and (|1.16|) we also see that the 
values of permanental polynomials close to the points i\i and — iji are 
non-trivially correlated in the limit of large N, as long as fi G (2,2). 
The latter fact could probably be related to the property of the per- 
manental roots to come in complex conjugate pairs (as is evident from 
coefficients of the permanental polynomials of Hermitian matrices being 
real). All this facts taken together make it rather plausible to expect 
permanental roots of large GUE matrices to accumulate asymptotically 
in the vicinity of the interval [—2i, 2i] of the imaginary axis. Below we 
provide an argument in favour of the validity of the following 
Conjecture: The normalized asymptotic limiting density of the per- 
manental roots of large GUE matrices in the complex plane z = x + iy 
is vanishing outside the interval [—2i, 2i] of the imaginary axis. Inside 
that interval it is given by the Wigner semicircular law. 

Our argument goes as follows. For z = x + iy define the function 
(1.17) *(x,y)= lim h\n\p(z)\ 2 ) GUE 

N—>oo iy 



According to the standard potential theory [3d] the expectation value 
of the limiting density p(x, y) of permanental roots in the complex 
plane is given by (p(x,y)) G uE = ^rA$(x, y), with A standing for the 
Laplacian in variables x and y (understood in the sense of distribu- 
tions). Unfortunately, performing the ensemble average in (|1.17J) is an 
outstanding problem. Assume however that the operations of taking 
the ensemble average and taking the logarithm commute in the limit 
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iV — > oo, i.e. 

(1-18) lim hln\p(z)\ 2 ) GUE = lim 1 In (\p(z)\ 2 ) GUE 

Such commutativity is indeed known to occur, for example, in the case 
of characteristic polynomials of GUE matrices[Zj, and one may hope 
it to hold for permanental polynomials as well. It would also imply 
the self-averaging property: the asymptotic root density is given by 
its expectation value. Proving relation (jl.l8|) for GUE permanental 
polynomials remains an outstanding problem. If however such a rela- 
tion holds, the function $>(x, y) can be found from (jl.l8j) . (jl.l5j) and 
(jl.lOj) by exploiting the known Plancherel-Rotach asymptotics for the 
Hermite polynomials 7Ttv(aO, see e -S- Eq.(7.93) of [18 . A simple calcu- 
lation yields for x ^ 

(1.19) $(x,y) = ^(y -ix) + ^(y + ix) 

(1.20) V(q) = -(q- vY - 4) 2 - In (q - vY - 4) + const 

8 

For x = Imq ^ the function ^f(q) is obviously analytic in the plane 
of complex variable q = y + ix, hence the function <&(x, y) is har- 
monic resulting in vanishing root density p(x, y) = 0. Thus, the only 
non-vanishing root density is possible for x — 0, i.e. along the imag- 
inary axis in the original complex plane z = x + iy. To calculate 
this density from <fr(x,y) given in ()1.19|) one again follows the stan- 
dard potential theory and evaluates the jump of the normal derivative: 
y)\ x ^o+ — ~§^®{ x i y)\x-*o- across the line x — 0. It immediately 
yields the Wigner semicircular law in the interval [— 2i, 2i] of the imag- 
inary axis, in agreement with the proposed conjecture. Similar, but 
slightly more technically involved calculations suggest validity of the 
semicircular density profile at the interval [— 2v2 i, 2y/2i] of the imagi- 
nary axis for the permanental roots of Gaussian Orthogonal Ensemble 
of real symmetric random matrices. 

The relation ()1.18j) is known to hold also for characteristic poly- 
nomials of Ginibre Ensemble and Circular Unitary Ensemble (CUE) 
matrices, see [27| and the discussion at the end of Sections |S] and El . 
Assuming its validity for permanental polynomials of those ensembles 
one can conjecture the ensuing asymptotic densities of the permanen- 
tal roots for those ensembles inside the unit circle x 2 + y 2 < 1 in the 
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complex plane. They are given by: 

(1.21) P(x,y) = - I , CUE 

7T (1 + x z + y 1 ) 1 

(1.22) p{ x iV) — j Ginibre Ensemble 

7T 

and vanishing density outside the unit circle. Whereas for Ginibre en- 
semble the conjectured density of permanental roots (|1.22|) is the same 
as that for the characteristic roots, the density profile (|1.21|) is obvi- 
ously radically different from the CUE eigenvalue density concentrated 
on the unit circle. Apart from investigating the validity of the above 
conjectures, any result on correlation patterns typical for permanental 
roots of large random matrices of various types are presently lacking. 
Finally, the question whether an accumulation of permanental roots in 
the vicinity of the imaginary axis is generic for more general ensembles 
of Hermitian random matrices is very intriguing and clearly deserve to 
be a subject of a separate study. 
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2. Permanents as contour integrals 

Let z = (zi, . . . ,Zn) t stand for iV— component complex (column) 
vector, and similarly define a (row) vector £* = (£ 1; . . . , £ N ), with the 
bar standing for the complex conjugation, and the star for Hermitian 
conjugation. Further denote z®£* the N x N matrix with entries z£j, 

i,j = 1,...,N, and let (£*z) stand for the scalar product Y^i^ii 
and £*Fz for the bilinear form Fij^Zj. 
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Lemma 2.1. The permanent of an arbitrary N x N matrix F = {F^} 
can be expressed through the contour integral as 

(2.1) PerF 



(2tt) 2JV 

jv 



| 21 |=1 J\ ZN \=lJ\^\=l J\£ N \=1 



x eX pTr[Fz®r]n^ 

k=i Zk £fc 

Proof. Obviously, Tr [F z ® £*] = £*Fz for any matrix F. Denote q = 
Fz, so that exp (f *Fz) = Ylf=i expfift. Use ^- §^ =1 1| exp £g = -q to 
see that the right-hand side of ()2.1j) is given by 

N N N , 

\^ ) J\zi\=l J\z N \=l i=1 j=1 k=1 Z k 



^ ji=l,..,iiv=l - / l^l= 1 "Vl =1 fc=l 2fe 

The last integral is obviously non- vanishing only as long as ji 7^ j% 7^ 
... 7^ j N , and the required relation immediately follows. □ 

Our way of representing permanents as contour integrals is of central 
importance for the rest of the paper. Despite its apparent simplicity, 
we were not able to trace such a formula in the available literature 
on permanents. In the Appendix we demonstrate how the formula 
(12. lj) generates a known integral representation of the permanent of a 
positive definite matrix as a multivariate Gaussian integral used earlier 
in the physical applications 54 . 

In what follows we will systematically use the short-hand notation 



(2.3) <t [...]V N (z,C 



N A N AC 



(2tt) 2jV J\ Zl \ =l J\z N \=i %l=i At N \=i k= i z l k=1 £ 
In particular, we will make use of the following 



k 



Lemma 2.2. Let f(z) be an entire function of the complex variable z 
represented by its convergent power series f(z) = Y^n=ofn zH - Then: 



(2.4) I N {f)= (b /(£*z) T) N (z, £*) — /, 



' A' 
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Proof. Use z k = e i9k , £ k = e iqpk , k = 1, . . . , N, where < 9 k , ip k < 2tt, 
so that 

i N ( f ) = (f\ 2N r r e -*™-*>f (j2e^A n^n^ 

\ 71 / Jo Jo \ k ) k k 

Shifting k ^ a k = 9 k — ipk, integrating out ipk and further reinstating 
the unimodular complex variable z = e iak we see that 

\**J J\ Zl \=l J\z N \=l J Z l Z N 

The integrations can be easily performed using the power series expan- 
sion, and the result follows. □ 

3. Mean value of the permanental polynomials for 
unitary invariant ensembles of hermitian matrices 

The main goal of the present section is to verify the formula (jl.lip . 
We need the following 

Lemma 3.1. Let V be any N x N matrix such that it has precisely one 
non-zero eigenvalue 7, the rest N — 1 eigenvalues being equal to zero. 
Then the following identity holds: 

(3.1) 7(/3;A,r) = / exp/5Tr [UAU*T] dji N {U) 

JU(N) 



(3-2) = E ( N N +n % ^rh n (Xi,...,X N ) 



n=0 

N 



m*- 1 ti iw a * - a,) 

where 

(3.4) h n (\i, . . . , Xn) = 2_s Xi x \i 2 . . . \i n 

l<i 1 <i 2 ...<i n <N 

are the complete symmetric functions. 

Remark: The integral in the right hand side of (j3.1)) is the well- 
known HarishChandra-Itzykson-Zuber (HCIZ) integral [313 E2|- I n the 
case of a general Hermitian matrix T with all distinct non-zero real 
eigenvalues 71, . . . , 7jv the result was found by Itzykson and Zuber 
in the form 

(3.5) /(ftA ,r) = r ~(n«!)^(f)- 
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In our case it is most convenient to exploit that fact that one can 
recover HCIZ expression by the Schur function expansion method 03] 
(also known as the character expansion, see [§]). In particular, the 
method does not require matrices A and T to be Hermitian, and we 
will indeed consider them to be general complex. 

Proof. Let us shortly remind the definition and basic properties of the 
Schur functions (see the book jJUj for more detail). A partition r = 
(ni, n 2 , . . . , um) is a finite sequence of integers, called "parts" , satisfying 
(jii > ri2 > ■ ■ ■ > tin), which is characterized by its weight |r| = Tij 
and its length l(r) equal to the number of non-zero parts. With any 
partition of the length l(r) < N one associates the Schur function 



dot {x^ +N - jXh 



(3.6) s r (x 1 ,x 2 , ■ ■ ■ ,x N ) 



\ - ■ " N 
det ( x- •' 



i 



which is a symmetric polynomial in x±, x 2 , ■ ■ ■ , x^ homogeneous of de- 
gree |r|. By convention, s r (xi, . . . ,xn) = for l(r) > N. 

For the simplest case of partitions of length one r = (n) (i.e. n\ = 
n, n 2 = . . . = tin = 0) the Schur functions coincide with the complete 
symmetric functions: 

(3.7) s (n) (xi, . . .,x N ) = h n (xx, . . .,x N ) 

defined in ()3.4|) . Using definition ()3.6|) . expanding the determinant in 
the numerator in entries of the first column and remembering 

det ( x f r ~"') = — x j) 



l<3 



we easily get a useful identity: 
(3.8) 



N ™n+AT-l 

^ Uj^i( x j ~ x i) ~ \ h n {xi,...,x N ), n = 0,1,2, 



. ir f_£^L = (0 1 n = _ (iV _i),_(iV_2 



For any N x N matrix M G Gl(N) one defines Schur functions of 
the matrix argument, s r (M), as 

(3.9) s r (M) = s r (mi, . . . , mjv) 

where are eigenvalues of the matrix M. The Schur functions of 
matrix argument are the characters of irreducible representations of the 
general linear group (and its unitary subgroup), and as a consequence 
have the important property of orthogonality. If T\ and r 2 are two 
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partitions, and Mi,M2 are two N x N matrices, then (see p. 445 of 

HOI) 

s ri (MiM*) 



(3.10) / s ri (VM 1 )s r2 (VM 2 ) d^ N (U) = 5 ri , r2 , 

JU(N) S ri {L Nj 

and 

s r (Mi)s r (M 2 ) 



(3.11) / 8 r {VM{irM 2 ) dfi N (U) (1 , 

JU(N) S r {LN) 

Using all these facts, one can calculate the integral in (|3.1j) by ex- 
panding the integrand in a series with respect to Schur functions, and 
performing the integration over the unitary group using ([3. II)) . After 
working out the coefficients of expansion the result is expressed as a 
sum over partitions [HJ EI] 

(3.12) I(/?;A,r) = J2(3 ni+ - +nN f[ J^~' )! -,, *r(AK(r) 

^ ■ L 4- (N + rii - i)\ 

for general matrices A and V . In our particular case 72 = . . . = 7at = 0, 
and therefore s r (T) = s r (7i) = 7™ for r = (n), i.e for partitions of the 
length Z(r) = 1, and s r (T) = for l(r) > 1. Using ([3.7)1 the expression 
([3.12)1 immediately takes the form of ([3.2)1 . and the equivalence to ()3.3|) 
follows after exploiting the identity ([3.8)1 . □ 

Now we are in position to prove Theorem ( II. lj) . 

Proof. Our starting point is the integral representation ([2.1)1 according 
to which the permanental polynomial p (/x) = Per^il^ — H) is given 

by 

(3,3) p M = /e X pMrz)e X p-r r[ Hz«r]^(z,r) 

According to ([1.6)1 , ([1.7)1 the ensemble average of this polynomial then 
amounts to evaluating the integral 
(3.14) 



{p(ji))=C N I e- NTrV ^A 2 (A) f[d\ t fe^J N (z,C)V N (z,C) 

J-co i=1 J 

where 

(3.15) J N (z,C)= [ exp-Tr[VAU*z<g>C]dfj, N (U) 

Ju(N) 
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For performing the latter integration over the unitary group we can 
use the Lemma ()3.1|) due to the fact that the (non-Hermitian) ma- 
trix z (g) £* has precisely one non-trivial eigenvalue equal to the scalar 
product (£*z). We find 



(3.i6) j n (z, a = jr ,S*^L (-ine*) n hn{Xi, . . . , a*) 

Substituting this back to 1)3.14)1 we at the next step should evaluate 
the contour integral according to (J2.4)) . For doing this we need to find 
N— th term in the Taylor expansion of 6^* z \Jn(z,^) in powers of (£*z). 
This is obviously given by 

(N-l)\^ 1 ^ (2AT-1)! f l)l+N _^ N _ lxl + N _ x 



1)! 1 V 

(2N - 1)! ^ n, ¥1 (A,- - Ai) i= _^_ i} (N -l)\(N + l- 1)! 

where we exploited the identity (|3.8jh which allowed us to replace the 
lower limit I = in the second sum in (|3.17|) with I = —(N — 1). Now 
the sum over I gives simply (// — Aj) 2iV_1 , and we arrive at the result of 
the contour integration in the form 

(JV-l)! A (p-X, 



(3,8) fe^J NM V^C) = g^jg^ _ Aj) 

Next step is to substitute ()3.18|) to (|3.14|) . and to use the symmetry 
of the integrand with respect to a permutation of variables Ai, . . . , Ajy. 
This gives 



(3.19) = CW / e-^^^-Ai)^" 1 ^! 

N N 



X 



/ e -^L a WJJ (Al _ Aj) JJ (A,-A,) 2 n^ 



j=2 2<j^i<N i=2 



where the particular value of the proportionality factor Cat is not 
needed for our purposes, and can be established whenever necessary 
( see below). 
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According to the classical theory of orthogonal polynomials |48[ IT%] 
122] . the integral in the last line of 1)3.19]) is proportional (see ()1.8|) ) to 
the monic orthogonal polynomial: 7Tjv_i(Ai) = A^ _1 + lower degrees, 
defined in (jl.9|) . The overall proportionality constant can be fixed 
by the obvious condition (p(/i)) = fJ, N + lower order terms, and we 
immediately arrive at the statement of the Theorem □ 

Corollary 3.2. For the Gaussian Unitary Ensemble (GUE) corre- 
sponding to the choice V(x) = x 2 /2 the relation 111. 11]) takes the form 



(3.20) (p(fi)) GU E = i 7Tiv(-i/i), 

Proof. The orthogonal polynomials are obviously Hermite polynomi- 
als. The corresponding monic polynomials TTk(x) have the following 
convenient integral representation, see e.g p. 53 of |25j : 

(3.21) 7r k (x) = (-i) k \[^e Nx2 l 2 r q k e~ N ^ +iNx «dq 

V 2tt y_ 00 
Substitute ()3.21|) to p. 11]) and change /i — A — > x, obtaining 



N— 1 / ju r°° r x 



Qn V 27T 



oo J — oo 



1)N -I* ^ x 2 ^ 1 -^ ( e~ N ^ N ^dq) dx 



axN*- 1 V 2vr J _^ dx N 



' i\iV roc jT. 



a N N N ~ 1 dx N 



The same method can be used to find the normalization constant from 
Ijrm and (I3~2T1) : 

(3 ' 23) aN ~W (N-l)\ ViV 

Comparing (13.2213.231) with (l3~2Tl) immediately verifies (]3~2U1) . □ 

In the next section we will be able to reproduce the above result 
for GUE case by a different method, which is essentially originating 
from applications of random Gaussian matrices in physics of disordered 
systems, see [23] and references therein. In this way we will be able 
also to get some insight into higher order correlation functions of GUE 
permanental polynomials. 
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4. Correlation functions for permanental polynomials of 

gue matrices 

The joint probability density of matrix entries for GUE N x N Her- 
mitian matrices H is given by 
(4.1) 



A' 



A? 



N 



V GUE (H)dH 



where Cque = 2 2 
stant. 



Cque Y\ ' 

i=l 

(N/ir) 1 



N 

, IV ^" n< v,/ ;/ n y// n 

i=l i<j i=l i<j 

A 2 is the corresponding normalization con- 



— TrA 2 



Lemma 4.1. The identity 

(4.2) {e- T ' HA )cu E 
holds for a general N x N matrix A with complex entries. 

Proof. The integral in the left hand side factorizes into a product of 
simple Gaussian integrals of the sorts 



g 2 a 



dHi, 



2tt 




iv exp 



—A 2 

2N " 



and 



e -NHij 11.. - 11,. Aji— "<v 1 >///:'• 



ij 



71 

N 



exp 



N 



A- A- 

l 3 3 1 



The product of all those factors together with the normalization 
constant Cque yields the right-hand side of ()4.2|) . □ 

The use of this lemma allows one to represent the ensemble averaging 
of the GUE permanental polynomial p(n) from ()3.13j) in the form 



(4.3) 



<PM> 



GUE 



exp 



Mrz) + — (rz) 2 



where we used that Tr(z ® ^*) 2 = (£*z) 2 . 

Let us now note that the identity (|4.2j) for the particular case of 
matrices H of the size lxl can be written in the form: 



(4.4) 



exp 



2N 




N 
2~H 



e 2 



N 9 

~ q - qa dq 



valid for any complex a and N > (we used an obvious replacement 
Ha —>■ q, An —>■ a and a simple rescaling of the integration variable 
with parameter N). We are going to use this formula for a = (£*z), 
in order to replace the factor exp [^(C 2 ) 2 ] m the integrand of ()4.3|) 
with the corresponding Gaussian integral. Such a procedure (and its 
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generalizations, see below (|4.12jl ) is customarily referred to as "the 
Hubbard- Stratonovich transformation" in physics [21]. 

Changing the order of integration, and exploiting (|2.1|) for the case 
F = (fjL — q)l N , i.e. 



j exp [(// - q)(Cz)} V(z, £*) = Per[(^ - q)l N ] = (ji-q) 



N 



we find 



(4.5) (p(»))gue=J^ I" e-^i 2 (fi - q) N dq 



oo 



which in view of (13.21)) is indeed equivalent to (|3.2(J|) . 

A nice feature of the present method is that it provides, after a sim- 
ple generalization, an access to higher-order correlation functions of 
permanental polynomials for GUE matrices. In fact, we will be able 
to reduce for any integer n = 1,2,... the general n— point correla- 
tion function of permanental polynomials of GUE matrices to the form 
([1.13)1 . and in this way to prove the Theorem (11.21) . Note, that the 
volume element <iq in (jl.13)) is defined in the way analogous to the 
measure dH in (|4.1|) . 

Proof. We consider in detail the calculation of the two-point correlation 
function {n = 2) defined according to ()2.1|) as 

(4.6) (p(/ii) p(fi 2 ))GUE = j j exp [//l(ClZl) + /i 2 (^2 z 2)] 

x (exp -Tr [H {z 1 <g> £ + z 2 <g> &)]) GUE V(z u £) V(z 2 , C 2 ) 

The case of general n can be verified precisely along the same route 
without any essential modifications. 

The ensemble averaging in the right-hand side of ()4.6)) is performed 
according to (|4.2j) . yielding 
(4.7) 

(exp -Tr [H ( Zl <g> Ci + z 2 ® &)]) GVE = exp —Tr ( Zl ® £ + z 2 ® £ 2 *; 
Now we notice that 

(4.8) Tr ( Zl ® £ + z 2 ® £*) 2 = (^ Zl ) 2 + (£ 2 *z 2 ) 2 + 2(£z a )(£z 1 ) 
and introducing 2x2 matrix 



,2 



(4.9) 



£l z 2 
£ 2 * Z 1 £ 2 * Z 2 



ON PERMANENTAL POLYNOMIALS OF CERTAIN RANDOM MATRICES 19 



we see that 

(4.10) (exp -Tr [H ( Zl <g> Ci + z 2 ® Q]) GUE = exp 



2N 



Traf 



Thus 
(4.11) 

{pMpM)gue 



j j e Mi(«zi)+^(e 2 *Z2)+^vTra 2 jy^ ^ jy^ 



At the next step we again exploit the Hubbard-Stratonovich trans- 
formation (cf. ()4.4j0 : 



(4.12) 



exp 



2N' 



-Tra 2 



e -fTrq 2 -Tr(qa) rf| 



where q is 2 x 2 Hermitian matrix 
(4.13) 



q 



9n 

(R) ■ CO 

912 -^12 



(R) , ■ (i) 

922 



q 



^ dgn rfg 2 2 ^91^ rfg^ 



The relation ()4.12|) follows immediately from ()4.2|) specified for 2x2 
matrices, after an appropriate change of notations. Substituting ()4.12|) 
back to (|4.11|) and changing the order of integrations, we arrive at 
(4.14) 

(pMpM)gue= J e-^ 2 dq j j e M^i)+M^ 2 )-Tr( q a) v ^ ^ p(z 

Now we note that 

Mf* z i) + ^2(^2) ~ Tr(qa) 
(4.15) = (/ii - gu)(£zi) - gi 2 (^ z i) - 9h(£i z 2) + (>2 - 922X^2) 



(//i-9h)1jv -912!^ 

—9121^ (^2 — 922) ljV 



Zl 

z 2 



which allows us to evaluate the contour integrals in ()4.14|) with help of 
the formula (|2.1j) . The result is simply 

(4.16) j j e W(^l)+M2(£ 2 *z 2 )-Tr(qa) jy^ ^ jy^ Q 



Per 



(A*i-9n)liV -912!^ 

— 9l2ljV (/^2 — 922)1^ 



Using that (i) the permanent of a block-diagonal matrix is equal to the 
product of the permanents of the diagonal blocks and (ii) the matrix 
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transposition does not change its permanent we finally arrive at the 
required identity: 

(4.17) tpM P M) G UE = I e~^ 2 [Per (M — q)] N dq 

where M = diag([i\, fi 2 ) ■ D 



According to the Theorem (|1.2|) the evaluation of the correlation 
functions of GUE permanental polynomials reduces to evaluation of the 
integral in the right-hand side of (jl,13|) . In the earlier studied case of de- 
terminantal polynomials performing a similar integration was achieved 
by shifting the integration variables M — q — > q, and passing to the 
eigenvalue decomposition q = U n qdU*, qd = diag(qi, . . . ,q n ), U n G 
U(n). As det q = detqd (and is therefore independent of U n ) one 
can evaluate the unitary-group integral with the help of HCIZ iden- 
tity which results in a great simplification. Unfortunately, permanents 
lack the nice invariance properties of the determinants, and for gen- 
eral n > 2 the exact evaluation of the corresponding integrals is still 
outstanding. The case n = 2 turns out however to be specific, and 
a simple modification of the procedure described above does yield the 
full solution of the problem. 

To this end, we note that for 2x2 Hermitian matrices q specified in 

Per (M - q) = - qn)(fi 2 - q 22 ) + M12 



;-l)det 



-111 + 9ll) -012 
-912 (A*2 - 922) 



Substituting this expression to the right-hand side of (|4.17jl ). and chang- 
ing qu — > —qu (which obviously leaves Trq 2 invariant), we therefore 
see that 



(4.18) (p(ft)pW) 



GUE 



-1 



e 2 



Trq 2 



det M - q 



N 



dq 



where M = diag(— /ii, ^2)- Now, as explained above, we can shift the 
integration variables M — q — >• q, and pass to the eigenvalue decom- 
position q = U2qdU2, qd = diag(qi, q 2 ), U2 G U(2). The measure is 
changed as dq — ► (rV 2 /47r)(q , i — q 2 ) 2 dqi dq 2 dfi 2 (U), and 



Trq 2 ^ Tr(M - q)^ = (^ + f4) - 2Tr MU 2 q d U* 2 + (q( + q z 2 
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Combining all these facts, we find that 
(4.19) (p(/ii) pM W = 



/oo /*oo 
/ e-^ { ^ + ^ 2) (gig 2 ) 7V (gi - q2) 2 dq 1 dq 2 I HCIZ (q 1 , q 2 ) 
■oo J —oo 



where we defined 



/ffc/z(?i,? 2 ) = Jexp-NTr (MU 2 q d U;) d/i 2 ([/) 

(4 ' 20) = Nfa + det {e N ^ e-*™ 

according to (J3.5j) as the standard HCIZ integral. After straightfor- 
ward algebraic manipulations the formulae ()4.19j) - (j4.2()j) ) give rise to 
the following final expression for the two-point correlation function: 
(4.21) 

{■nfn W,» 1 Hp+ ( TTiv(-^l) 

(pMP^W " (Z^) _ (-) det ^ ^ +1 (^ 2 ) 



where we have used the integral representation 1)3.21)1 for the monic 
Hermite polynomials. 

In particular, last expression proves a curious relation (J1.15)) between 
correlation functions of permanental and characteristic polynomials of 
GUE matrices, valid in view of the corresponding Brezin-Hikami result 

CQDJ. 

5. Two-point correlation function of the permanental 
polynomials for matrices from the circular unitary 

Ensemble 

Consider the space U (identified with the unitary group U(N)) of Nx 
N unitary matrices, satisfying U _1 = U*, with the probability measure 
on U being chosen according to the Haar's measure on the unitary 
group. This condition defines the Circular Unitary Ensemble (CUE). 
With the permanental polynomial defined as p(n) = Per(fil^ — U), it 
is immediately clear that any non-trivial correlation function must have 
the form of a product of an even number of permanental polynomials, 
precisely half of the factors involving matrices U*, that is 



=1 / CUE 

In what follows we will address the simplest non-trivial case of the 
two-point function ( n = 1), which according to (|2.1j) has the following 
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integral representation: 

(5.2) x ( e -2>[u(-A®G)+u-(«B®a)]^ V(z A ,e A )V{z B , 

To perform the ensemble average of the exponential factor in the inte- 
grand above we make use of the following 

Theorem 5.1. For two general complex square N x N matrices A and 
B define 

(5.3) F(AB*)= [ expTr [AU + U*B*] dfj, N (U) 

JU(N) 

Suppose AB* has rank 1. Then for N > 2 holds 

(5.4) F(AB*) = (N - 1) / (1 - t) N ~ 2 I (2\/itf) eft 

Jo 

where v 2 is the only non-zero eigenvalue of AB* and 

°° z 2j 

is the modified Bessel function. 

Proof. The theorem is just a special rank— 1 case of a more general 
statement proved for any rank 2m < N in |2Zj using the Schur function 
expansion. It is intimately related to the so-called "bosonic colour- 
flavour" transformation introduced by Zirnbauer in jSS]. For rank m = 
N the function F(AB*) was earlier calculated in see also j^j. □ 

According to (|5.2j) we are interested in a particular case A = -za <8> 
C A and B* = -z B ® so that AB * = (G z b) z A ® wni ch is 
obviously of rank 1, with v 2 = (C a z b)(Cb z a)- Using (|5.4j) and the 
expansion (J5.5|) it is convenient to perform integration over t explicitly, 
and find: 

(za\ / P -Tr[u(»A®tt)+u«(«B®ft)]\ _/ v iy [(G z b)(G z a)] J 
1 ; \ /crus 1 j\(N-l+j)\ 

Substituting this expression to ()5.2|) we need to evaluate the integral 



(5.7) 

= jj e V iA) ^^'M [(Ca^)(Cb^a)} j V( Za ,Ca)^b,Cb) 
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This is most easily done by first considering the generating function 



(5.8) I0 A \^), Px , P2 ) 

= j j e [^(az A )+^T(^Z B )+ P l(a Z B)+ P2 (^ZA)] p ( ^ £* ) p (zs; £* ) 

from which ()5.7j) is obtained by a simple differentiation, and then taking 
limit pi = P2 = 0. Noticing 



(5.9) //< A >(Gz A ) (&z B ) +Pi(Gz B ) +p 2 (£W) 

/i (A) ljV piljv \ f Z A 



B) 



p 2 l N Ai (B) liv / V z ^ 
we can apply the Lemma ()2.1|) to ()5.8|) and find 



(5.10) I^ A \^), Pl ,p 2 ) 



Per 



/i (A) pi x 1 ' 

P2 



where we have used that the permanent of a block-diagonal matrix is 
equal to the product of the permanents of the diagonal blocks. Now a 
simple differentiation yields the required result: 



(5.11) h(u^ A \d B )) = I (N-j)i^ » ) > 

V ; JVP P ; [ 0, j >N 

Combining all these results, and changing (N — j) — > j for convenience 
we arrive at the final expression for the two-point correlation function 
Q: 

N i TB)\j 

(5.12) ( P{I ^)W^) CUE = - D! g flfa I ,■), 

A direct check shows that the above expression can also be written in 
the form of an integral: 

(5.13) (p(^ A) )p(i^)) ={N-1) I (l-t) N - 2 {ii^J^)+t) N dt 

\ I CUE 





This latter form is especially suitable for extracting the asymptotic 
behaviour in the limit of large N by the Laplace method. In particular, 
for the integral is dominated by the vicinity of the point 

of maximum of the integrand t = t m = [1 — | /x | 2 ) /2 for < 1, and 
by the lower limit of integration t = for |/i| > 1. A straightforward 
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calculation then shows that 
(5.14) 

W) - Ijm iln<WM)l'> OT - { ^f^M 1 

Under assumption that the relation (|1.18|) holds for CUE matrices, the 
expression (|5.14j) for <3> p yields after applying to it the Laplace operator 
a two-dimensional density profile in the complex plane fi = x + iy given 
in the equation (jl.21j) . 

Since this result is very different from the density of (unimodular) 
eigenvalues of the CUE matrices, let us briefly show how to recover the 
latter density in a similar approach. The analogue of ()5.13|) was found 
in |2Zj, and reads 
(5-15) 

(d(^>) W^)) CUE = (n + i) J™ (1 + l t)N+ M A) ^ + tf dt 

where d(fi) = det (//Ijv — U) is the corresponding characteristic poly- 
nomial. In contrast to the previously considered case, in the limit of 
large N the integrand for li^ = li^ = ll does not have a maximum 
in the domain of integration, and thus the integral is dominated by the 
vicinity of the upper limit t = oo for < 1, and by the vicinity of the 
lower limit t — for \fi\ > 1. This yields 

(5.16) MW 2 ) = M ^(WP)\% UE = { 'fJ.jVi 

We see, that the function y) is harmonic everywhere in the com- 
plex plane \i = x + iy, except the unit circle \fi\ = 1. This means that 
the density of eigenvalues of CUE matrices is indeed vanishing every- 
where except \z\ = 1, in agreement with obvious unimodularity of the 
corresponding eigenvalues. The discontinuity of the normal derivative 
of the "potential" $i across the unit circle yields the constant density 
on the circle, as expected. 

6. Correlation functions for permanental polynomials of 

GlNIBRE MATRICES 

The joint probability density of matrix entries for complex Ginibre 
N x N matrices Z is given by 

N N 

(6.1) V GUE (Z)dZ dZ* = C Gin H e~ N ^ H dZf 3 dZ\ 3 

where Coin — (N/7r) N is the corresponding normalization constant. 
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Lemma 6.1. The identity 



(6.2) 



-Tr(ZA+BZ*) 



Gin 



e N 



Tr AB 



holds for any two general N x N matrices A, B with complex entries. 

Proof. The integral in the left hand side factorizes into a product of 
simple Gaussian integrals 



e N 



The product of all those factors together with the normalization con- 
stant Cque yields the right-hand side of ()6.2|) . 



□ 

Again, as in the case of CUE the nontrivial correlation functions for 
Ginibre ensemble must contain even number of factors, with half of 
them being permanental polynomials of Z, the rest being permanental 
polynomials of Z*. More precisely, we will prove the following 

Theorem 6.2. For any integer n = 1,2,... the 2n— point correla- 
tion function of permanental polynomials of Ginibre matrices Z has 
the following integral representation: 
(6.3) 

Q* 



\k=l 



p /4 



Gin 



-NTrQ 2 



Per 



M A 

Q 



M 



B 



dQ dQ* 



(A) 



and~M B = diag(n[ B \ ; 



In this expression = diag(fi\ 
are n x n diagonal matrices and Q is a general n x n complex matrix, 
with the volume element dQ dQ* defined in the way analogous to the 
measure dZdZ* in \6. 

Proof. We consider in detail the calculation of the two-point correlation 
function (n = 1) defined according to (j2.1|) as 

(6.4) (p(A* (A) )pP 5 5)) ojB = j je^[^ A \e A ZA)+Jl {B \eB^)] 
x (exp -Tr [Z z A ® & + Z*z B ® C B ]) Gin V(z A , &) V(z B , Cb) 

The ensemble averaging in the right-hand side of (J4.6|) is performed 
according to ()6.2|) . and using Tr (z A ® £4) (zb ® £ B ) = (£aZ B ) (£s z a) 
yielding 
(6.5) 

(exp -Tr [Z z A ® L + z * z b ® Cb] ) Gin = exp 



N 



(£>b) (&z a ) 
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Introducing the notations 

(6.6) P AB = (ftz B ) , Pba = (? B z A ) 

we notice that the relevant form of the Hubbard-Stratonovich transfor- 
mation in this case is given by (cf. 1)6.2)1 ) 

(6.7) e MPABP BA ] = E J J e -NQQ + QP AB+ QP BAdQ (R) dQ( I) 

where the integral over Q is in the plane of a general complex variable, 
Q = Q R + iQ 1 . Substituting this back to ()6.5|) and then to ()6.4|) and 
changing the order of integrations, we have 

(p (v {A) ) p¥^)) Gin = 7 / / e ~ mq d Q d Q 

x I ^exp [/i (A) (az A ) +7i (B) & B ) + Q(£>b) + Q(€bZa)] V(z b ,C b )V(z b ,Cb) 
Further noticing that 

/i (A) (e>A)+/I (B) (GzB)+Q&B)+g(^Z A ) = ( ^qi n N ^b)^ ) ( II 

the contour integral can be performed using (j2.1|) . yielding finally 
P (^))H^)) = -/ fe- N ^Per( ^ Jj" ) dQ dQ 

I Gin 71 J J \ l^i-N H y ' i-N J 

Per 



(6.8) = ~ f f e~ NQQ ' X1 A 



7T 

in accordance with ()6.3|) . 



dQ dQ 



The proof for general n follows precisely the same lines with minimal 
changes. The place of two parameters Pab and P B a defined in (|6.6|) is 
taken by two n x n matrices Pab an d Pba with entries 

(6-9) [Pab] 13 = (&j*B d ) , [PbAh = (Cb^a,) 

where i, j = 1, ...,n, and the corresponding Hubbard-Stratonovich 
transformation is 

(6.10) e ^Tr[P AS P flA ] = J e -iVTr[Q*Q]+Tr[QP AB +Q*P SA ] d Q rf Q* 

□ 
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Note: The problem of evaluating the integral in the right-hand side 
of (|6.3|) for general n > 1 is still outstanding. For n = 1 the calculation 
is elementary: 

(v^ {A) )JW^)) Gm = "J J e~ N ® Q [/i (A) /J ( ' B) + QQ] N dQ dQ 

poo 

(6.11) = N / e""* + R) dR 

Jo 

k=0 

The integral form (|6.11|) is most convenient for extracting the large-N 
asymptotics by the Laplace method. For /i^ = fj,( B > = fi the integrand 
has a sharp maximum around R = 1 — |/i| 2 as long as < 1 and is 
dominated by the lower limit R = for > 1. This gives: 

(6.13) * Gi „(^| 2 ) = ^Ihn iln{b (/J )|^ Gi „ = { 1^-^ W^ 1 

Under further assumption of validity of (|1.18|) such an expression im- 
plies the uniform density of permanental roots inside the unit circle, as 
conjectured in (|1.22|) . 



7. Correlation functions for permanental polynomials of 

GOE matrices 

The joint probability density of matrix entries for GOE N x N real 
symmetric matrices H is given by 

N N 

(7.1) v GO E(a)dH = cgoe n e_f hi n e ~ NHfj n dH » n m ^ 

i=l i<j i=l i<j 

where Cgoe = 2~^~ (N /t[) n ^ n+1 ^ a is the corresponding normalization 
constant. 

It is straightforward to verify the following 
Lemma 7.1. The identity 

(7.2) ( e ~ TrUA ) GO E = eM TrA2 + TrAAT ] 

holds for any general N x N matrix A with complex entries. 
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The use of this lemma allows one to represent the ensemble averaging 
of the GOE permanental polynomial p(fi) from (J3.13j) in the form 
(7.3) 

N \ 2 „ N N 



(p(h))goe = j> exp 



m 



-rl 



where we used that Tr(z <g> £*)(z ® £*) T = J2i £i Si z f- 

Following the same strategy as for GUE case, we use the appro- 
priate gaussian integrals ("Hubbard-Stratonovich transformation") to 
linearize the quartic terms in the exponential. The first term is identi- 
cal (up to the factor 1/2) to that in the GUE case, and we can use (|4.4|) . 
The second term in the exponential is dealt with exploiting (|6.7|) . with 

the correspondence Pab ~^ Si£i> Pba — ► ^2i z f- Exchanging the or- 
der of integrations, and remembering that P(z, £*) is essentially given 

^ WF? Uti if nf=i f (see Q), we arrive at 




(7.4) (pQi))aoE =\/^ I dq- / e^^dq.dq, 



oo 



X 



27T) 2 ^| = l^| = l Z 2 f 2 



-I iV 



Noticing, that the result of the contour integration is independent of 
(as any analytic factors depending on z 2 or £ could be replaced 
by zero's term in their Taylor expansion) we see that second line in 
(|7.4jl obviously produces (fi — ^q) N ' , and after simple manipulations 
and exploitation of (j3.21|) we arrive at the following 

Proposition 7.2. The expectation value of the permanental polynomial 
for the Gaussian Orthogonal Ensemble (GOE) is given in terms of the 
monic Hermite polynomial \3.21)) by 

i N 

(7.5) (p (^)}goe = 2jv/5 TfvHv^/i) = i N (d {-i^))goe , 

where again d(p) = det (/iljv — H). The second equality follows 
from comparison with the equation (48) of |T2"] . 

The same method can be used, mutatis mutandis, for evaluating 
higher correlation function, although actual calculations are much more 
cumbersome than for the GUE case. Below we restrict ourselves to 
considering explicitly the two-point correlation function, n = 2. As the 
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starting expression we use the analogue of ()4.6j) : 

(7.6) <pMpM)goe = j> j> exp [/ii^zx) + /x 2 (^z 2 )] 

x (exp -Tr [H ( Zl ® £ + z 2 ® £ 2 *)]> G o£ P ( z i> £) P ( z 2> 8) 

where the ensemble average is now performed using the identity ([7.2)1 . 
with the role of A played by the matrix z x ® £* + z 2 <8> £ 2 . 

According to (|4.8jl . for such a matrix TrA 2 = Tra 2 , where a was 
introduced in (|4.9j) . We also notice that 

TV TV TV TV TV TV 

Tr AA r = X) S< X) 4 + E & E 4 + 2 X ?m?« E 

i=l i=l i=l i=l i=l i=l 

Combining Hubbard- Stratonovich transformations (|4.12|) with (three 
times) (|4.4|) we see that the ensemble average featuring in ()7.6|) can be 
represented as 

e -Tr[H( Zl ®£r+z 2 ®£ 2 *)]\ = f e -TVTrq 2 -Tr(qa) rfq f e ~N{Wqi +9292+23393) 

/ gob 7r 3 y y 

xe lELi [ft Ef =1 Ef =1 Ef =1 ^«+WEf=i ? M ?u] dq 1 dq 1 dq 2 dg 2 dq 3 dq 3 

Substituting this back to 1)7.60 . changing the order of integration, and 
suppressing the terms which do not contribute to the value of the con- 
tour integral, we obtain after a simple manipulation: 

(7.7) (p(jJLl)p(jJb))GOE 

2N J e-^+a**) [J(q,g 3 ,^i,/i 2 )f rfq^: 



71 

where Trq 2 = q\ x + g 22 + 2q l2 qu and 

(7.8^(^93,^3;^,^) = -^/ 4 4 4 e 5 3 ? 1 f 2 +« 3ZlZ2 



3 



(27t) 4 J| Z1 | =1 J| {i [ =1 J|„[ =1 J| & [=1 



x exp 



/"i - 011 -012 \ / «i 

-012 A*2 - 022 / V ^2 



j fc=i fc fc=i 4fc 



= [(^1 - 01l)(^2 - 022) + 0l 2 012 + 0303] 

Changing finally gn — > — 0n we arrive at the final expression 
(7.9) (p(ft) P« W = ^ / e-^ +a *») 



X [(AH + 0il)G"2 - 022) + 012012 + 0303]^ dqdq 3 dq 



3 
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Comparing this result with the equation (43) of [12] we infer the validity 
of the relation: 

(7.10) (p(A*i) pM)goe = (d{Wi) d{-i^ 2 ))GOE , 

which is precisely the same as for GUE case, (jl.l5|) . 



Appendix A. Permanents as multivariate Gaussian 

INTEGRALS 

We start with verifying the following 

Lemma A.l. Let F > be positive definite N x N Hermitian matrix, 
a = (di, . . . , ajv) T , b = 6/v) T , s = (si, . . . , sn) t all be N— com- 

ponent complex vectors, and ds*ds = Y\^ =1 dRe{si)dIm{s-i). Then the 
following identity holds 

(A,) e— * = detp/ e — <^W* 

Proof. Representing F = UfU* in terms of the diagonal matrix of the 
eigenvalues f = diag(/i, . . . , /jv) > and the unitary matrix U of the 
eigenvectors, and introducing new variables of integration q = U*s so 
that ds*ds = dq*dq, the AT— fold gaussian integral in the right-hand 
side decomposes into the product of A" simple Gaussian integrals of 
the type flB77|) : 

N N 

CLi hi 



(A.2) fl I e-^'-^-^'dqidqi = f[ ef 

i=l J i=l 



where we introduced the notations a = U*a and b = U*b. We further 
find that 



N 1 - f N 1 \ 

(A.3) J2j~ a ^ = E U ^7 U ^ 

i=l Ji jk \i=l Ji J 



a*F~ 1 b 



which proves the statement. □ 

Proposition A.2. The permanent of a positive definite Hermitian ma- 
trix F can be represented in terms of a multivariate Gaussian integral 
as 

1 r N 

(A.4) Per F = — - / e - s * F " s J] (s fSi ) ds*ds 

J i=i 
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Proof. According to the contour representation of the permanent (J2.1j) 

(A.5) PerF = j e^ Fz V N (z,C). 

We now replace the exponential factor in the integrand with its repre- 
sentation according to the Lemma (jA.ljl : 

(A.6) e ? * Fz = det F 1 J e^ F ' ls -^^ds*ds 



N 

(A.7) j e-^-^V N (z,C)=U 



Substituting this back to (|A.5|) . changing the order of integrations and 
straightforwardly performing the contour integration according to 

N 

[Si*, ; 

i=i 

we immediately arrive to (53) . □ 

Note: The formula (jA.4|) is in fact a variant of the Wick theorem 
for the Gaussian integrals. 

Further note that any positive definite Hermitian matrix F can be 
represented as F = EE*, where E = Uf 1 / 2 in terms of the eigenvectors 
and the corresponding eigenvalues of F. This allows us to arrive at the 
following result, which seems first to appear in |54j : 

Corollary A. 3. The permanent of a positive definite Hermitian matrix 
F can be represented as 

N 

(A.8) PerF = / e^ JJ ([Ev*] i [Ev].) dw*dw 

J i=i 

Proof. Introduce in (|A.4|) a new integration variable v according to s = 
Ev. The integration volume is changed as ds*ds = det E det F*dv*dv = 
det Fdv*dv, and the statement follows. □ 
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